Dressed-state approach to a simple radiating atom 
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We consider a system consisting of an atom in the dipole approximation, coupled to the elec- 
tromagnetic field. Using recently introduced renormalized coordinates and dressed states, we give 
a non-perturbative solution to the atom radiation process. From this formalism a non-exponential 
decay law naturally emerges. However, at least for low frequencies and weak coupling, there is only 
a very small, practically non-detectable, deviation from a purely exponential decay, in agreement 
with all previous results in the literature. 
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I. INTRODUCTION 

Actually, apart from computer calculations, the most currently used method to solve coupled equations describing 
the physics of interacting systems is perturbation theory. This is a common feature to different branches of physical 
sciences, such as celestial mechanics, field theory and statistical physics. The method, originally introduced to deal 
with orbital problems in celestial mechanics, has been extremely successful since then, giving remarkably accurate 
results in Quantum Electrodynamics and in Weak interactions. In high energy physics, asymptotic freedom allows to 
apply Quantum Chromodynamics in its perturbative form and very important results have been obtained in this way 
in the last decades. 

However, in spite of its wide applicability, there are situations where the use of perturbation theory is not possible, 
as in the low energy domain of Quantum Chromodynamics, where confinement of quarks and gluons takes place. 
In this branch of physics, few analytical approaches in the context of Quantum field theory are available up to the 
present moment (in this particular subject the interested reader can find a recent account in [1]). There are also 
other situations, in the domain of Cavity Electrodynamics and Quantum Optics, where perturbation methods are 
of little usefulness, for instance, the observation of resonant effects associated to the coupling of atoms with strong 
radio- frequency fields. The theoretical understanding of these effects on perturbative grounds requires the calculation 
of very high-order terms in perturbative series, which makes the standard Feynman diagrams technique practically 
unreliable [2]. Methods to non-perturbatively deal with such kind of systems, have been introduced originally in Refs. 
[3] and [4] , and since then they have been used to investigate several situations involving the interaction of atoms and 
electromagnetic fields [5-7] . 

In this paper we give a non-perturbative treatment to the system composed of an atom (approximated by an 
harmonic oscillator) interacting with the electromagnetic field. In the dipole approximation, we reduce the problem 
to a system of a mechanical oscillator (the bare atom) coupled to the infinite set of the harmonic modes of the 
(bare) field. The emergence of an exponential decay of an excited state of the oscillator and also renormalization 
aspects of this problem have been thoroughly discussed in the literature. See for instance Refs. [8-15]. Here we use 
an alternative approach, based on precisely defined renormalized coordinates introduced and employed recently to 
describe the coupling of a mechanical oscillator with a scalar field in previous papers [16-18]. In fact, in the quoted 
papers these coordinates have been denoted as dressed coordinates but, as explained in the paragraphs following Eq. 
(26), we presently think that the denomination of renormalized coordinates is more appropriate for them. It should 
be noticed that our definition of dressed states is not the same of the collective dressed states introduced in previous 
works in the literature [3,4]. In particular our dressed states are not eigenfunctions of the interacting Hamiltonian. 

To describe the radiation process, having as initial condition that only the mechanical oscillator (the atom) be 
excited, we do not proceed as usual, considering the interaction term in the Hamiltonian written in terms of the 
bare oscillator coordinates and field modes as a perturbation, which induces transitions among the eigenstates of 
the free Hamiltonian. We introduce renormalized coordinates for, respectively the mechanical oscillator and the 
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electromagnetic field modes. In terms of these new coordinates, dressed states can be consistently defined in section 
III and in terms of these a non-perturbative treatment of the radiation process is possible. 
Unless explicitly stated all along the paper we use natural units c = l,Ji= 1. 



II. THE COUPLED SYSTEM 



In this section we review the basic aspects of the method of reducing the system composed of a harmonic oscillator 
(the atom) interacting with the electromagnetic field to a system of coupled harmonic oscillators. Readers familiar 
with this derivation can skip the calculations and go directly to the comments at the end of this section. 

In the Coulomb gauge, 

= 0, V.A = 0, (1) 
the Lagrangean for the oscillator-electromagnetic field system is given by 

,2 



(2) 



where the dots mean derivative with respect to time, 

J = eS{x — q)q (3) 

and e is the electric charge of the oscillator. We consider the particle-field system in a cube of side a, that is, 
—a/2 < {x,y,z) < a/2 and impose periodic boundary conditions for the field: A{—a/2,y,z,t) = A{a/2,y, z,t), 
A{x,—a/2,z,t) = A{x,a/2, z,t) and A{x,y,—a/2,t) = A{x,y,a/2,t). At the end we will take the continuum limit 
a — > 00. 

Expanding the electromagnetic field in real modes we have 

A{x, t) = J2^k {x)qk (i) , (4) 

k 

where the Ufc(x)'s form a real orthonormal basis and satisfy 

-V^Ukix) = Lol . (5) 

Replacing Eq. (4) in Eq. (2), using the identity V. = G.V X F — F.W X G and the orthonormality relations 

of the Ufe(^)'s, we find easily, 

m ■_2 vauP' _2 1 

' k 

For all dynamical purposes the above Lagrangian is equivalent to 



_rru^^2 _ _Y^(^^^2 _ ^2^^ 2^ ^ ^j^Uk {q)'q.Qk ■ (6) 



L = - - ^ Z] (^fc ^ - ^l^k ^) - eY^Uk{q)q.qk > (7) 



k k 



since the difference between Eqs. (6) and (7) is a total time derivative. 
Solving Eq. (5) with periodic boundary conditions we find 



where V = , Wk = \k\ and 



/2 -* / 2 ^ 

- sin(fc.x) , w^^^ (f ) = _ cos(fc.x) , (8) 



271" 

k = — {l,m,n) , Z, m, n = ±1, ±2, ±... (9) 



Notice that there are positive and negative values of I, m, n. However only one half of the total set of possible values of 
A;'s will be admitted since sm{k.x) and cos{k.x) are not, respectively, linearly independent of sin{—k.x) and cos(— fc.a?). 
Replacing Eq. (8) in Eq. (4) and using the transversal condition V.A = we find 

£gfe(«)(t) = 0, a = 1,2, (10) 

where qk^^^{t) and qk'''^\t) are associated, respectively to and u''^\x). Prom Eq. (10) we can write 

Qk^^^t) = q^^^ki + Qk2^k2, where Cki and efc2 are the unit orthogonal polarization vectors. Next we consider 
the dipole approximation, that is, wc assume that the electromagnetic field docs not changes significantly for points 
around the oscillator equilibrium position g = 0. This allows us to make in the last term of Eq. (7) the substitution 
Uk{q} ~ Wfc(O). Prom Eq. (8) we find Ufc^'(O) = 0, from which we conclude that, in the dipole approximation, the 
modes qk^^^{t), do not interact with the oscillator. On the other hand uf'\o) = Also, if we choose one of 

the polarization vectors to lye on the plane defined by k and q we can see that only the mode corresponding to this 
polarization vector interacts with the oscillator since the other mode will be orthogonal to q. If we choose e^i as 
orthogonal to q, then the terms different from zero in the last sum of Eq. (7) will contain terms proportional to 
q-^2kq^h ~ sin^/j^gfe'? where 6k is the angle between q and k. Then Eq. (7) can be cast in the form, 

^=^1 ^1 - 2 W ~ '^fe^fc) ~^iv l^^^^^^qkq , (11) 



where we have written only the field modes that interact with the the oscillator (without labels, since there is just only 
one) and we have suppressed the vectorial symbol in the oscillator coordinates since now only appears its modulus. 
Prom Eq. (11) we obtain the Hamiltonian, 

k k k 

Notice that in the above Hamiltonian the last term proportional to q^ is infinity. This infinite term guarantees the 
positivity of the Hamiltonian. In the case of an oscillator coupled to a massless scalar field, considered in Refs. [16-18], 
this term is absent at the beginning and is introduced afterwards by renormalizing the oscillator frequency. Indeed, 
also in the present situation we can define a bare frequency lob, such that, 

a;|=a;2 + — 2_]sin6ife. (13) 



mV , 

k 



By performing a canonical transformation, 



p-*\/rnpo and qk^—, Pk ^ -^^kqk (14) 



in Eq. (12) we get 



where 



H=^i + ^^l + \Y.ipl+ ^l<ll) - E ^'^^^ ' (15) 

k k 



Ck = ey' ^ smOkL^k ■ (16) 

The Hamiltonian given by Eq. (15) is almost identical to the corresponding Hamiltonian in the scalar field case [16-18], 
with the differences that the physical oscillator frequency arises naturally (that is, we do not need to renormalizc the 
frequency oscillator), the coupling involves an angle 6k and the sum is in three dimensions rather than in one. We 
have to mention that a Hamiltonian of the form given in Eq. (15) has been used many times in the literature as a 
simplified model to the quantmn Brownian motion [19,12] and also to study the decohcrencc problem [10]. However 
we would like to remark that in all early treatments no use of renormalized coordinates (to be defined below) has 
been done in the way we employ them. Basically we introduce our renormalized coordinates in order to guarantee 
the stability of the ground state of the oscillator (the atom) in the absence of field quanta (the vacuum), (see Refs. 
[16,17] for details). If these coordinates were not introduced, as explained below, the oscillator in its ground state 
would decay to other states, in contradiction with experimental facts. 



III. RENORMALIZED COORDINATES AND DRESSED STATES 



We can diagonalize the Hamiltonian given by Eq. (15) as in Refs. [16-18] (see also appendix of Ref. [20]), defining 
= J2r ^liQr, Pn = J2r ^li^r, M = (0) k), where {V^} is an orthogonal matrix. We obtain, 



(17) 



The matrix elements are given by 



1 + V ^ 



(18) 



and the normal frequencies Qr, are given as solutions of the equation 



Using Eq. (13) and (16) we can write Eq. (20) as 



(19) 



(20) 



(21) 



In order to define the renormahzed coordinates and dressed states we start from the eigenstates of our system, 
represented by the normalized eigenfunctions, 



(Q)=n 



ro((3) 



where stands for the n^-th Hermite polynomial and 



(22) 



(23) 



is the ground state eigenfunction of the system. As in Refs. [16-18] we introduce renormalized coordinates q'^, 
/x= 0, fc = l,2,..., 



(24) 



where uj^ = {u!.{u!k}). In terms of the renormalized coordinates we can define our dressed states as the complete 
orthonormal set of functions [16], 



ro(g') 



where ro(g') has the same functional form as ro{Q), that is, 

ro(g')oce~^SM"''< 



(25) 



(26) 



It is worthwhile to notice that our renormalized coordinates are new objects, different from both the bare coordinates 
q and the normal coordinates Q, and have been used for the first time in Refs. [16-18]. In particular, the renormalized 
coordinates and dressed states although being collective objects, should not be confused with the normal coordinates 



Q, and the cigenstates (22). Differently from the eigenstates (22) which are stable, our dressed states are all unstable, 
except the ground dressed state obtained by setting {m^ = 0} in Eq. (25). 

The idea, introduced in Ref. [16], (for a recent clear explanation see Ref. [21]) is that the dressed states given by Eq. 
(25) are supposed to be the physically meaningful states, instead of the ones written in terms of the bare coordinates 

that appear in the original Lagrangean. This can be seen as an analogous of the wave function renormalization in 
quantum field theory, which justifies the denomination of renormalized to the new coordinates q' . Thus, the dressed 
state given by Eq. (25) describes the atom in its mo-th excited level and ruk field quanta (photons) of frequencies 
(jjk ■ It should be noticed that the introduction of the renormalized coordinates guarantees the stability of the dressed 
vacuum state, since by definition it is identical to ground state of the system. The fact that the definition given by 
Eq. (24) assures this requirement, can be easily seen replacing Eq. (24) in Eq. (26). We obtain TQ{q') oc ro(Q), 
which shows that the dressed vacuum state given by Eq. (26) is the same ground state of the interacting Hamiltonian 
given by Eq. (17). 

The necessity of introducing renormalized coordinates can be understood by considering what would happen if we 
write Eq. (25) in terms of the bare coordinates 5^. In the absence of interaction, the bare states are stable since they 
are eigenf unctions of the Hamiltonian obtained from Eq. (15) by setting Ck = 0. But when we consider the interacting 
term they become all unstable. With the excited states there is no problem since we known from experiment that 
this is indeed the case. On the other hand, we know from experiment that the atom in its ground state is stable, in 
contradiction with what our simplified model for the system describes in terms of the; bare; coordinates. Obviously 
our simplified model is only approximate, the correct theory that describes the system atom-electromagnetic field 
is Quantum Electrodynamics (QED), but a rigorous, non- perturbative description using QED would be extremely 
complicated. So, if we wish to have a non-perturbative approach in terms of our simplified model we have to modify 
something in order to remedy the above mentioned problem. The solution is just the introduction of the renormalized 
coordinates q'^, as the physically meaningful ones. On the other hand if, in terms of creation and annihilation operators, 
we use the so called rotating wave approximation [22] for the interacting term in Eq. (15) (the Lee-Friedrichs model), 
wc would remain, neglecting part of the interaction Hamiltonian, only with terms involving ajofe and a[,ao [this is 
not the case for the whole Hamiltonian (15) which we consider here]. In this situation there is no problem with the 
stability of the atom ground state described in terms of the bare coordinates g^, since in this case the action of the 
interacting term on the bare oscillator ground state (obtained from Eq. (25) setting q'^ — q^ and ~ 0), gives zero. 
That is, the free oscillator ground state still remain as the ground state of the interacting system. This is the reason 
why in early treatments (see for example [23] and references therein) of the simplified model presented here, there has 
been no need to introduce renormalized coordinates. 

We can write the Hamiltonian (15) in terms of dressed coordinates. Replacing Eq. (24) in Eq. (15) we get the 
renormalized Hamiltonian, 

^ = \ ("oo(Po)^ + Pooco^iqof) + i'^kkip'kf + Pkkcoliq'kf) + ^ {aokP'oP'k + PokOJo^^Jkq'oQk) 

k k 

+ {^krk2PkrPk2 + ^krk2i^k^(^k2<lkA2) ' 

ki>k2 

where 
and 

/3m. = -z^E^«- (29) 

Expression (27) seems complicated and without clear physical significance. In particular, it appears that the dressed 
atom and the dressed field modes are in the same foot of equality, that is, the dressed atom, described by q'^ and Pq, 
appears to play no central role as in the original Hamiltonian. However this is not the case. Sinc;e the expressions for 
tQ and t]., given by Eqs. (18)-(19), are very different, the renormalized coordinate and momentum p'^ still play a 
central role, because they are multiplied respectively by aoo and /3oo that contain terms of the type (tfj)^, see Eqs. 
(28)-(29). On the other hand the dressed coordinates and momenta pj. are multiplied by Ukk and I3kk than contain 
terms of the type tyff^. The interacting terms in Eq. (27) involve both dressed coordinates and momenta different 
from the original Hamiltonian, that involve only explicit interaction among the coordinates go and Qk- The last term 



(27) 



(28) 



in the first line of Eq. (27) expresses explicit interaction between the dressed atom and the dressed field modes. The 
second line in Eq. (27) expresses the explicit interaction among the dressed field modes. However the couplings of 
these two type of interactions arc very different, the dressed atom-photon and dressed photon-photon interactions 
contain, respectively, terms like tjjt^ and t^, that are very different. 

We can show that the Hamiltonian written in dressed coordinates, given by Eq. (27), approximately decouples into 
the dressed atom and the dressed field modes. That is, we show that in Eq. (27), a^^ ^ Oi^i,\^i^i, and /S^^ ^ fifiv\ix^v 
For this end note that for = (or equivalently for e = 0, where e is the electric charge) =9^- In this case, 
Q^MM = /^MA* = 1 ^■'^d "^i^Ut^i^ PiJ.u\iJ.^u = 0. Then, because e < 1, we must have a^^ = 1 + 0(e), and = 1 + 0{e) 
whereas a^i,\n=^v = 0{e) and I3^i,\^j^^i, = 0{e), showing that ^ ct^^^ and /3^^ ^ Pfi^v This shows that the 
introduction of the renormalized coordinates allows, approximately, to divide the system into the dressed atom and 
the dressed field modes. 

The decomposition given in Eq. (27) shows that the renormalized Hamiltonian takes into account explicitly some of 
the implicit interactions that are preset in the original Hamiltonian. For example, a photon-photon interaction term, 
of the type given by the second line of Eq. (27) can be obtained from the interaction term q^)qk^ that appears in the 
original Hamiltonian, at second order in perturbation theory. Also we have to note that in Eq. (27) the " uncoupled" 
terms, are not identical to the uncoupled ones in Eq. (15), the differences are just self-interacting terms. From these 
discussions we can see that our dressed states, defined by Eq. (25), are not free, they are collective states, that take 
into account the interactions. The dressed atom, for example is the bare atom "surrounded" by bare (virtual) photons. 
The same picture applies to the dressed photons. That these claims are correct can be seen naturally, writing any 
dressed state given by Eq. (25) in the basis of bare states, the bare states being identical to the states given by Eq. 
(25) with the dressed coordinates q'^ replaced by the bare ones q^. 

With the formalism of dressed coordinates and dressed states we can compute the probability that the atom 
(understood as the dressed mechanical oscillator q'^) prepared at f = in its first excited state, still remain at the 
time t in the same state, or from it the probability of decay of the atom from its first excited state. For this end we 
consider the initial dressed state V'ioo...(9')i which only the dressed atom is in the first excited level, all the others 
oscillators being in the ground state. Then it is shown in Ref. [16] the following expression for its time evolution, 

^ioo...((z',i) = ^/o^(t)^o...ovo...(a') (30) 

where 

k^{t)=Y.^itle-^''''- (31) 

Prom Eq. (30) we see that the initially excited dressed atom naturally distributes its energy among itself and all other 
dressed oscillators (the dressed field quanta) as time goes on, with well defined probability amplitudes given by the 
quantities fo^{t) in Eq. (30). From Eq. (30) we obtain for the probability amplitude that the dressed mechanical 
oscillator (the atom) remain in its first excited state the following expression, 

/oo(t) = E(*o)'e-'""* . (32) 

r 

As already mentioned above, we do not make explicit use of the concepts of interacting bare oscillator and field, 
described by the bare coordinates go and {qu}- We have introduced renormalized coordinates and {g^} for, 
respectively the dressed atom and the dressed field modes. In terms of these new coordinates dressed states can be 
rigorously defined. A non-perturbative approach to the radiation process is possible in terms of the dressed states 
for the system contained in a cavity of arbitrary size, including microcavities [17]. Also free space emission starting 
from periodic boundary conditions, which we examine in this paper can be treated using our dressed states. This 
unified treatment is not possible in the conventional framework. Our dressed states Eq. (25) are collective but non 
stable states, linear combinations of the eigcnstatcs (22) defined in terms of the normal modes. The coefficients of 
these combinations can be calculated and explicit formulas for these coefficients for the interesting states are given in 
Refs. [16,17]. Moreover, as noticed above, they have the very interesting property, used in Refs. [16-18] for the case 
of a sc;alar field, of distributing the energy initially in a particular dressed state, among itself and all other dressed 
states with precise and calculable time-dependent probability amplitudes. We choose these dressed states as physically 
meaningful to the study of the radiation process. 



IV. THE RADIATION PROCESS 



Now, to compute Eq. (32) in the continuum limit a — > oo, let us define the complex valued function, 



(33) 



Notice, from Eqs. (21) and (33), that the fi^'s are given by the real roots of W{z), and that in the continuum hmit 
these real roots will extend over all the positive real axis in the complex ^-plane. Deriving Eq. (33) with respect to z 
and comparing with Eq. (18) we find. 



2z 



W'{z) 



Replacing Eq. (34) in Eq. (32)and using Cauchy theorem we get 



1 

iri 



dz 



W{z) 



(34) 



(35) 



where C is a counterclockwise contour in the ^;-plane that encircles the real positive axis. The integral in Eq. (35) can 
be computed choosing a contour C that lies just below and above the real positive axis. Below the real positive axis 
we have z = fl — ie and above z = fl + ie. Then we have 



/( 



1 r 



Wirt-ie) W{n + ie) 

The next step is to compute W{^1 + ie) and — ie). Using Eq. (16) in Eq. (33) we get in the continuum hmit 



(36) 



W{z) = z'^-ij^ + 



= .2 — W + 



4m 



ir/2 



ddk 



d,^ 



37r2 



m 



dk 



(37) 



where we integrate over the 6k angle only up to n/2 instead to tt, since as we remarked above, only one half of 

the modes arc linearly independent. The integral over k in Eq. (37) is linearly divergent. Then in order to make 
meaningful the above equation we will need to regularize it. We do this by analytic regularization, that is, we introduce 
the integral 



7(a) = ^l' 



a-l / 

Jo 



dk 



fc2 



{k'^-z 



2\a 



(38) 



which is well defined for a > 2 ; then we can give a meaning to the integral in Eq. (37) by analytically extending 
the integral given by Eq. (38) to a = 1. The parameter /i in Eq. (38) is introduced in order to maintain the correct 
dimension of the original integral. Using Cauchy theorem we obtain for I {a), 



I{a) = iTrn'^-\a - 1) 



d° 



I{a) = -iTrfj,"-\a - ly. 



dw"'-^ {w + zY 



Im{z) > , 
Im{z) < . 



(39) 
(40) 



Then by analytic extension we find 7(1) = mzj^ for Im{z) > and 7(1) = —iTrz/2 for Im{z) < 0. Replacing these 
values in the integral of Eq. (37) we find 



w{n + ie) = _ 

W{il - ie) = _ 



Gttto 



+ 0{e), 
+ 0{eY 



(41) 
(42) 



Substituting the above expressions in Eq. (36) we get in the e — > 0"*" hmit, 



f00it)^-J^ ^» (^2 _ ^2)2 +,2^6 ^ (43) 

where we have introduced the time dimensional parameter r = e^/(67rm). The integral given by Eq. (43) can 
not be evaluated analytically. We can check the validity of Eq. (43) taking t = and noticing that for tuj ^ 1 the 
integrand is sharply peaked around Q = lu. Then we obtain from Eq. (43), 

/oo(0) ^ / dx^- — 

= 1 , (44) 

as expected from the definition given by Eq. (32). One can check the above result also for any value of r using Cauchy 
theorem for i = in Eq. (43). We see from Eq. (43) that the probability |/oo(i)P will not be a pure exponential. 
However, at least for low frequencies as we will see below, there will be only a small deviation from an exponential 
law. 
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FIG. 1. Plot on the vertical axis, of the probability that the oscillator remain excited at time t in units cot = s. 

Considering t > we can integrate Eq. (43) numerically and from the result compute the probability P{t) — |/oo(^)P 
that the atom still remain excited at time t. For instance, taking tlu — 0.05 and using for e and m respectively the 
proton charge and the mass of the hydrogen atom, we get uj « 5.9 x 10^°/s. For these values we obtain the plot 
showed in Fig. 1 where we can see an almost exponential behavior for the decaying process. 



V. CONCLUDING REMARKS 



In this paper we have analysed a simplified version of an atom-electromagnetic field system and we have tried to 
give the more exact and rigorous treatment we could to the problem. We have adopted a general physicist's point of 
view, in the sense that we have renounced to approach very closely to the real behavior of a complicated non-linear 
system. As a counterpart, an exact solution has been possible. We have used a formalism [renormalized coordinates 
and dressed states) that allows a non-perturbative approach to the time evolution of the system. In particular, we 
have applied this formalism to study non perturbatively the time evolution of the excited atom. In the free space 
limit, we get the result that the renormalized coordinate approach describing the atom modified by the presence of the 
field in an indissoluble way, gives an exact result for emission in free space, generalizing the well known exponential 
decay law. 
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